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ABSTRACT 

We applied the ATON code to the computation of detailed grids of standard (non- 
rotating) and rotating Pre-Main Sequence (PMS) models and computed their adia- 
batic oscillation spectra, with the aim of exploring the seismic properties of young 
stars. As, until now, only a few frequencies have been determined for ~ 40 PMS stars, 
the way of approaching the interpretation of the oscillations is not unique. We adopt 
a method similar to the matching mode method by Guenther and Brown making use, 
when necessary, also of our rotating evolutionary code to compute the models for 
PMS stars. The method is described by a preliminary application to the frequency 
spectrum of two PMS stars (85 and 278) in the young open cluster NGC 6530. For 
the Star 85 we confirm, with self-consistent rotating models, previous interpretation 
of the data, attributing three close frequencies to the mode n=4, 1=1 and m=0,+l,-l. 
For the Star 278 we find a different fit for the frequencies, corresponding to a model 
within the original error box of the star, and dispute the possibility that this star has 
a Teff much cooler that the red boundary of the radial instability strip. 

Key words: stars: oscillations - stars: evolution - stars: rotation ~ stars: variable: 
other - stars: pre-main-sequence 



1 INTRODUCTION 



The existence of pulsating stars among intermediate mass 
Herbig Ae/B e stars (1.5 < M/M© < 4) was originally 



suggested by iBregei] (llOTd) and confirm ed by the pulsa- 
tion of HR 5999 l|Kurtz fc Mulled Il999l ) and HD 104237 
l|Donati et al.|[l997l V Considering that young PMS stars like 
the Herbig Ae/Be stars evol ve through the instabilit y region 
of post-main sequence stars. iMarconi fc PaUal (|l998l ) studied 
the radial pulsating properties of young stars, defining their 
instability strip (IS) for radial pulsations. Following work 
focused also on the search of further candidates, and today 
several objects have bee n found and are discussed in the 
astronomical literature (|Marconi et al ] I2OO4I : iRipepi et al.l 
[2006 a: Zwintz 200i). Most of these stars have more than 
one pulsation frequency, they look like PMS counterparts of 
the main sequence (MS) or early post-MS variables of the 
5Scuti type. It is now established that the PMS variables 
can also pulsate in non radia l p-modes (|Ripepi et al. I l2006bl . 
I2OO7I : IZwintz. Guenther fc Weiss 200^). Unfortunately, not 
many frequencies are identified in observations made from 
the ground, and their analysis then lacks firm identification 
of modes. Nevertheless, new observations are beginning to 
be available and more are expected from the program of cur- 



rent (COROT) and up coming (KEPLER) space missions. 
To complicate the things these stars have inherited an angu- 
lar momentum from the previous protostellar phase, which 
explains why most of these stars rotate. Recently, the ef- 
forts to measure stellar rotation periods — or at least the 
projected rotational velocities v sin ^ — for many stars of dif- 
ferent type have been intensified (|Rover et al.ll2002al lbh. In 
particular projected rotation al velocities in very young clus- 
ter stars have been measured (He rbst et al ]|2005l :l iPetit et al.l 
I2002,) with the aim of providing constraints in the mod- 
elling of the evolution of angular momentum during the very 
early phases of stars formation and PMS stages. In partic- 
ular the study of rotational velocity in the Orion complex 
l| Wolff et al.ll2004D has shown that their variation with the 
evolutionary phase, from the convective PMS to the Main 
Sequence, is compatible only with decoupling of core and 
envelope rotation when the stars begin evolving along the 
radiative tracks. 

Since the stellar structure is significantly modified by rota- 
tion, rotation will affect not only the age and mass of the star 
derived from the location in the HR diagram, but also the 
determination of its physical parameters and of oscillation 
frequencies in particular. Till now the mode identification, 
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that is the determination of radial order and spherical har- 
monic of the observed oscillation spectra, has been based 
on non rotating evolutionary models for the computation of 
theoretical frequencies, in some cases only noticing that the 
observations present a split, probably due to the rotation. 
As more and more new data from ground and space based 
observations are likely to become available in the near fu- 
ture, we decided to begin developing theoretical tools for 
the interpretation of the PMS variability, by using standard 
and rotating PMS models computed wi th the code ATON 
l|Ventura et al.ll 19981 : iMendes et al.lll999l ). and deriving their 
seismic properties by rneans of the adiabatic oscillation code 
LOSC jScuflaire et all l200iD . In Section 2 we discuss our 
theoretical tools: the input physics of the models, and the 
frequency computation. The method adopted for comparing 
observed and theoretical frequencies is described in Section 
3. We test the method, in Section 4, on Star 85 of NGC 
6530, and compare our results with those of lGuenther et aD 
I2OO3) (hereafter G07). We find that rotating models are con- 
sistent with the hypothesis of the split in a triplet of the 
frequency at 180.26 /iHz. In Section 5, we also discuss the 
comparison with Star 278 of NGC 6530 whereas in Section 
6 we summarize the main findings of this paper. 



2 THEORETICAL TOOLS 
2.1 The evolutionary code 

To build the grids of models we have used the code ATON 
l| Ventura et al.lll998l) in its standard v ersion for asteroseis- 
mic applications l|D'Antona et al.ll20()5l ) and in its rotational 
version, including stellar rotation according to the formula- 
tionby [EndaUsSofia ( 1976)_ as described in.Mcndcs et aTl 
l|l999l ) and in lLandin et all feood). The program has be en re- 
cently updated I Ventura. D'Antona fc Mazzitellill2007^ ■ and 
we describe in the following its main physical and numerical 
inputs. 



2.1.1 Input micro-physics 

The r adiative opacities are ta ken from Iglesias fc RogersI 
l[l993), and extended by the iFerguson et al.l ta- 
bles in the low temperature (T^ 10000 K) regime; the 
conductive opacities, harmonically added to the radiative 
opacities, are taken from Poteckhin (2006, see WEB page 
www .io ffe. rssi.ru / astro/con duct/). The OPAL equation of 
state of lRogers et al.l l|l996l ). overwritten in the pre ssure ion- 
ization regime by the EOS bv lSaumon et all (Il995l ). is used. 
The EOS is extended to the high-den sity, high-temperature 
domai n according to the treatment bv lStolzmann fc Blockerl 
l|2000h . A detailed description of the EOS will be presented 
elsewhere (Ventura fc Mazzitelli, in preparation) 



2.1.S Convection treatment 

The convective regions can be described either by 
the traditional Mixing Length Theory (MLT) approach 
l|Bohm-Vitense I Il958l). or following t he Full Spectrum of 
Turbulence (FST) (jCanuto et al.lll996l ') prescription. All the 
models presented in this paper have been calculated us- 
ing the FST treatment. Since this work is focused on PMS 
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Figure 1. Theoretical HR diagram for the PMS evolution phase. 
Evolutionary tracks calculated by ATON for th e labelled masses 
arc reported. Overplotted is the IS computed bv lMarconi fc Fallal 
(il998|). 



evolution, no extra-mixing was considered, though the code 
presents the possibility of allowing some extra-mixing (both 
in the instantaneous and the diffusive modality) from any 
convective border. 



2.1.3 Chemical composition 

All th e models were computed with the iGrevesse fc Noels! 
(| 19931 ) solar mixture of heavy elements, and the mass frac- 
tion of (X,Y,Z) is (0.70,0.28, 0.02). The initial mass fraction 
of Deuterium is X(D) ~ 2 • 10"^ 

2.1.4 Rotation 

Rotation was modelled according to the treatment by 
lEndal fc Sofial (|l976l ). This approach accounts only for the 
hydrostatic effects of rotation, neglecting the internal angu- 
lar momentum redistribution. Three rotation schemes are 
currently implemented in the code, namely i) rigid rotation 
of the whole star, ii) local angular momentum conservation 
everywhere, iii) local angular momentum conservation in the 
radiative regions and rigid rotation of the convective zones. 
The initial total angular momentum of the star is provided 
as a physical input. In this wo rk we used the iii) option in 
agreement with the analyses of I Wolff et al.l (|2004l ) . The ini- 
tial angular momentum is chosen to reproduce the observed 
surface angular velocity of the star at a given position in the 
HR diagram. 

2.2 The model grids 

Fig. [T] shows a selected set of our evolutio nary tracks from 
1.6 to 4Mq and the radial IS found by iMarconi fc Pallal 
( 1998) for the first three radial modes. In detail it represents 
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the instability strip limits between the second overtone 
blue edge and the fundamental red edge. The red edges 
for non radial oscillations all lie at Tofi larger than this, 
so the fundamental red edge represents the lo wer effective 
temp erature for any kind of PMS oscillation (|Unno et alj 
1 19891 ). This limit, however, depends on the efficiency of 
convective treatment, being cooler for a smaller efficiency 
of convection, as the interaction with deep convection 
quenches pulsations. T he red edge of the radial IS of 
iMarconi fc Pallal (|l998l ) was obtained with 1/Hp = 1.5 that 
reproduces the solar radius in their models. Decreasing the 
mixing lengh parameter from 1.5 to 1, the red edge becomes 
cooler by 200K (Marconi & Palla, private comunication) . 
Models were computed from 1.5 to 4 M© with a step 
AM=0.02 M0. The evolution begins in the Hayashi track, 
and the calculation ends when the star has reached the 
Main-Sequence. We memorize the structural quantities 
needed to perform the computation of the oscillation 
spectrum. The step in effective temperature of this grid is 
A Toff ~ 50K. For rotating models we follow a different 
approach. As the moment of inertia of the model changes 
along the evolution, together with the total radius and 
the internal mass distribution also the angular and surface 
velocity change along the PMS track, increasing from the 
youngest model (at the lowest examined Tcff) to the oldest 
model (at the hottest Tofr). We build up a sub-grid of 
rotating models after having compared observed frequencies 
with those of the non rotating grid and having extracted 
an initial list of best fit models. 



splitting Si^n between consecutive m components within each 
(l,n) multiplet,i.e: 

fi,n,m = .fi,n.o + m/3i,„n/27r (1) 

where /3, the coefficient of the rotational splitting, is calcu- 
lated directly by LOSC for each model. We stress that, as 
a novelty in the present computation, fi^n.o is the eigenfre- 
quency of (l,n) mode of the stellar rotating model, structure 
of which depends on the rotation rate. 

In Fig. [2] we show how these frequencies change for a model 
of the same mass (=2.40Mq) and effective temperature ('^ 
7140K) but with three different values of equatorial velocity 
of rotation (v=0 (NULL), v ~ 50 km/s (INT) and v~ 80 
km/s (IfIGH)); the figure shows the percentage variations 
of frequencies with 1=0,1,2 and n^26. The general trend is 
that the differences increase when the rotation velocity in- 
creases. It is interesting to note that for n>3-4 the frequen- 
cies of the rotating stars are smaller than those computed 
for non-rotating models (<Q.2% in the INT case and <1% 
for HIGH rotation velocities). For small n, on the contrary, 
the oscillation frequencies of the rotating models are larger 
then the non rotating ones, up to 2-3%, although the per- 
centual differences globally decrease for increasing n. We see 
that there is some scatter in the points at low n (see, e.g., 
the scatter in dots and triangles at n = 6,7,8). This is part 
of the problem appearing in the computation of frequencies 
for small n in the LOSC code. The problem will be care- 
fully investigated in the near future, as PMS stars oscillate 
mainly at this low n modes. 



2.3 Adiabatic oscillation code 

We use the LOSC oscillation code bv lScuflaire et all (|2007l ) 
to calculate the frequencies of radial and non radial oscil- 
lation modes for each grid point . We limit the spherical 
degree to 1=1,2,3,4 and cover the interval in angular fre- 
quencies between 0.6 and 60 times the dynamical timsQ- AH 
these oscillations are computed with the standard surface 
boundary condition SP/P -I- (4 + Lj^)Sr/r = 0, where uj = 
27r f Tdyn is the dimensionless angular frequency. 
Beyond the radial order n and the degree I another number 
characterizes the stellar oscillation frequencies fn,i,m'- the az- 
imuthal order mG[-l, 1]. The rotation breaks the azimuthal 
symmetry and removes the 21-1-1 degeneracy in m. 
In our modelling, we consider the symmetrical part of the 
centrifugal distortion directly in the computation of rotating 
stellar models. Then we take into account the first term of 
perturbation. The perturbation method at the first ordeiQ in 
the rotational velocity SI predicts an equidistant frequency 



'^dyn = y Qjj 

^ In fact recently, iLovekin fc Deupree I | |2008l ) using detailed 2D 
stellar models and a 2D finite difference integration of the lin- 
earized pulsation equations to compute non radial oscillations 
have definitively demonstrated that the eigenfunction can be ac- 
curately modelled using the perturbation theory and a single 
spherical harmonic, if rotation is slow enough (v^q <90 km/s). For 
these velocities is justified to neglect perturbation orders greater 
than one. 



3 COMPARING OBSERVED AND 
PREDICTED FREQUENCIES: THE 
METHOD 

The methodology that we will adopt to model the sparse 
radial and non-radial oscillation spectra observed in pul- 
sat ing PMS stars was origin ally suggested and developed 
by ICuenther fc Brownl l|2004h and used recently by G07 to 
model the oscillations o f five pulsating star s in NGC 6530. 
As recently reviewed bv ICunha^ elaD|2003), there are sev- 
eral other methods for mode identification of the frequencies 
that constitute the observed spectra of a given star and/or 
to constrain the mass and age of this star. These meth- 
ods use at the same time radial and non radial frequencies. 
Among them, one of the most popular consists in the use of 
the asymptoti c theory and of the l arge and small separation 
in particular (|Ruoppo et al.ll2007l ). This is a good strategy 
if the star oscillates in the asymptotic regime, or when the 
oscillation data are not well determined but, for the higher 
quality data expected from satellites, e.g. from COROT, a 
more flexible choice is to use all the diagnostic informations 
obtained fro m the data, and in partic ular all the discovered 
frequencies (jCuenther fc Brownll2004l ). 

Indeed, direct fitting of the frequencies is a relatively simple 
way to find a good first solution to the problem of mode 
identification, and this will become an ever and ever better 
tool with higher quality data. In principle, one could pro- 
ceed with seismic modelling comparing all values for (Z, n) 
to each of the detected frequencies. The value of degree I is 
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Figure 2. Left panel: PMS tracks corresponding to M=2.4 Mq and computed with different rotation velocities (see text). Right panels: 
differences (in percentage) for each pair of n and l(0-asterisk, 1-open dots, 2— filled dots) between frequencies of rotating models and 
NULL ones for the models labelled with an asterisk in left panel. 



usually limited (1 ^ 4jf| and m=0 is assumed when no obvi- 
ous evidence of rotational splitting is present in the observed 
spectrum. Given the A'^ frequencies of a particular star, into 
our grid of the oscillation spectra we look for models having 
^1, where: 



1 (Uob 



\2 

^mod.i ) 



N ^ 

i = l 



obs.i 



(2) 



t'obs,! and Vmod,i are respectively the observed frequency and 
the corresponding model frequency for the ith frequency, 



(Toftj j and Cmod i the associated uncertaintiefl In general 
the search for the best-fit model (which has the spectrum 
that minimizes the x'^ ) is done only in a sub-grid made of 
models within the 3cr error box in Toff and L of the star. 
If the program does not find any models with x^^li that 
means that not all frequencies are matched within the er- 
rors, it discards one by one the non-matched frequencies, 
computes again the x^ for the remaining frequencies and 
searches again for models with x^^l- The models with the 
lower value of the x^ is the best-fit model. 
When possible, we use the observed frequencies according 
to the order with which they are extracted from the Fourier 



^ This consideration is based on a geometric statement. In fact 
in distant stars, only the large scale structures can be detected 
and the sensitivity of brightness variations is restricted to modes 
up to 3 or 4 node lines at the surface. 



* For the available observations the model uncertainties are neg- 
ligible, as they are an order of magnitude smaller than the ob- 
servational ones (see Guenher & Brown, 2004 for a discussion on 
this topic) 
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Figure 3. The position of Star 85 in tlie HR diagram (see text). 
Solid lines are evolutionary tracks calculated with the ATON code 
with no rotation. Dashed lines represent tracks (with mass be- 
tween 2.00 and 2.50 and pass in mass of 0.02) calculated 
with option of rotation with an initial angular momentum which 
gives models with v=10km/s near the triangles. 

analysis of the data, and we try to fit the frequencies follow- 
ing this order. The order assigns to the first frequency the 
highest probability to be a real one, because it is the first 
detected by the Fourier analysis and, as a consequence, it 
deserves the largest number of options in the space of pa- 
rameters. The same reasoning applies to the following fre- 
quencies. Instead, the theoretical frequencies for each model 
are examined in order of increasing I (and for each I of in- 
creasing n) so that the possible matching begins by compar- 
ison with radial oscillations first. This approach gives then 
a higher probability to have a matched frequency with 1=0. 
Afterwords, the search of the closest model frequency for the 
second one in the list (which, being second has a lower am- 
plitude of oscillation) is now restricted to model frequencies 
that have not been already selected, and so on. 
If there is an evidence for the rotation of the star, for ex- 
ample an evident triplet in the observed spectrum, or other 
rotational indicators from photometry or spectroscopy, we 
calculate a sub-grid of rotating models in a box around the 
best fitting model(s) previously obtained. The rotational ve- 
locity chosen is that consistent with the mean separation 
observed, with the use of Eq. 1. Then, in the same way, we 
look again for the model(s) that minimize through the 
models of the new sub-grid. 

In the next section we apply this method on two test stars. 



4 EXAMPLE-1: STAR 85 OF NGC 6530 

Among several oscillati ng PMS stars fo u nd in the young 
cluster NGC 6530 by IZwintz fc Weis3 (|2006l ) (hereafter 



Table 1. Observed frequencies of Star 85. Those labelled with (I) 
are found by ZW06 while (II) marks the new frequencies found 
during the recheck performed by G07; components of the split are 
labelled with an asterisk. In third, fourth and fifth columns are 
reported 1, n and m of the best— fit modol(see text). 



No. 


Frequency (/iHz) 


1 


n 


m 


fl(I) 


180.26* 


1 


4 





f2(I) 


147.02 





4 





f3(I) 


179.65* 


1 


4 


-1 


f4(I) 


122.58 





3 





f5(I) 


360.48 





12 





f6(II) 


180.95* 


1 


4 


1 


f7(II) 


133.38 


2 


1 





f8(II) 


167.47 


2 


3 





f9(II) 


327.31 


2 


10 






ZW06), Guenther and collaborators have studied in detail 
the observed spectrum of five stars and were able to con- 
strain the models using their mode matching method for 
several stars. In particular, for Star 85 (WEBDA 53) several 
frequencies are available, and they show an evident triplet 
probably due to rotation. For this star, ZW06 identified five 
frequencies (labelled with I in Table 1). Noticing that two of 
them are too much close together, G07 have rechecked the 
data analysis using more recent reduction algorithms. In this 
way they found the lacking frequency of triplet (which ap- 
pears as a 1 day~^ alias) and three additional frequencies 
(all labelled in Table 1 with II). Using ±0.5 /iH^ as ob- 
servational uncertainty and keeping only the frequencies of 
the first run (without f3 which belongs to the triplet), G07 
find several models with ^ 1 , but only one within their 
observational box around the position of Star 85 in HR dia- 
gram and chose this one as the best fit model. Three of the 
frequencies were fitted with 1=0 modes and those suspected 
to be splitted were fitted with an 1=1 p mode. Among the 
frequencies emerged during the recheck they preserved their 
best fit model and found that the new frequencies lie very 
close to 1=2 p modes, except for the highest one (see Fig. |3] 
and[Sl). 

4.1 First step: non rotating models 

We have performed a search similar to that by G07, trying 
to fit the observed frequencies, in the order for which they 
are given in Table fl We begin with those found by ZW06, 
i.e. with those labelled with I in Table 1. Assuming for lumi- 
nosity and effective temperature of Star 85 the mean values 
given in ZW06 we consider an error of ±0.05 in logTcff and 
±0.3 in logL/L© (see Fig. |3]|. Inside this box we search for 

^ This error correspond to one over the time lengh of the obser- 
vations. 

^ We have chosen a different order with respect to Table 2 of G07; 
in fact for the oscillation modes of the first run we have preserved 
the original order, so the amplitude of oscillations decreases going 
from the first to the last frequency as explained in the previous 
section. Unfortunately this is possible only for the modes found 
in the first run, for which amplitudes and detail of the Fourier 
transformations are given and not for the other four frequencies 
which are only provided in increasing order in G07. 
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BEST-FIT MODEL 

M = 2.26M0 

T^f,(non rot) = 6843K 
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Figure 4. Comparison between tlie observed spectrum of Star 85 and the spectrum of the best fit model (filled triangle in Fig. |3j. Solid, 
dotted, dashed and long dashed lines connect theoretical frequencies with respectively 1=0, 1, 2, 3, 4. In the small box we report a zoom 
around the frequency 180.26 /jHz in order to underline the split do to the rotation of the star. 



our best-fit models. This box is larger than those used by 
G07 as we have noticed that their best-fit model lies just 
on the cool edge of the observational box. Only three mod- 
els fit within the errors all fouiQ frequencies (x^ =S! !)• In 
Fig. |3] we report these models with triangles; the filled one 
is that with the lowest value of x^y our best fit model 
which is practically overlapped to the model of G07 (filled 
circle). Also the mode identification is the same, with the 
three frequencies f2, f4, f5 matched with (l,n)=(0,4), (0,3), 
(0,12), and the splitted one (fl) with (l,n) = (l,4) (see Fig. 

Alternatively, by fixing these four pairs of (l,n) to the ob- 
We have discarded 13 which is a splittet frequency. 



served frequencies, the models with the lowest (including 
the three ones previously found) all lie on a diagonal of HR 
diagram (see Fig. |3]) as found for the case B by G07. This 
can be understood in the following way: by choosing these 
pairs (l,n) we assign a large separation (= — /„_i^;), 
which in asymptotic regime d epends linearly on lu minosity 
and temperature as shown by IRuoppo et al.l l|2007h . 
Concluding, using similar hypotheses, we obtain the same 
result of G07, i.e. the same best-fit model. The same mod- 
els are selected, when we search among all the grid models. 
Adding the frequencies found during the recheck of the data 
doesn't change the model with lowest x^ but in this case 
the lowest x^ is not ^1. Looking in detail to the compari- 
son, the main reason for this is that the three new frequen- 
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Figure 5. Echelle diagram for star 85, using the tlieoretical fre- 
quencies of tlie best fit model shown in Fig. |4] The parameters 
used to reproduce this plot are also shown. 

cies (f7, f8 and f9 in particular) are not perfectly matched 
(A/ > aobs from the nearest ones with 1=2). This problem 
was also found in the G07 investigation. A possible interpre- 
tation is given in the next subsection. 
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Figure 6. The position of Star 278 in the HR, diagram (see text). 
Solid lines are evolutionary tracks calculated with the ATON code 
with no rotation. 



4.2 The second step: rotating models 

Now we implement the comparison with rotating evolution- 
ary models. Using the observed splitting Af=0.643±0.09 
fiHz {= ), Eq. 1 and the radius of the best-fit model 
found in the previous section we derive a rotation rate of 
Veq ~ lOkm/s. We calculated a sub-grid of rotating models 
choosing an initial angular momentum that allows to obtain 
the predicted velocity near the position in HR diagram of 
the best-fit model previously found. 

The model grid is restricted to 2^ ^2.5 and has the 
same step in Tetr and mass of the non-rotational grid. The 
rotational velocity is Veq ~ 10 km/s around log Teff=3.83 . 
We compared again observed and calculated spectra looking 
for a good match in terms of ■ Such a slow rotation veloc- 
ity modifies very little both the structure of the star and its 
position in HR diagram, so we obtain identical results. Al- 
though we find the same three best fit models by considering 
the fit of the first frequencies only, when we add the three ad- 
ditional frequencies the rotating models match them better 
but continues to lie la away from the nearest theoretical fre- 
quency. Maybe it would be useful to recheck the data of the 
second round. An alte rnative explanation for this d i screp - 
ancy is the following. lOaszvnska-Daszkiewicz et al.l (|2002l ) 
tested the results of non adiabatic models on the P Cep stars. 
One of the effects of an even slow rotation, is to influence 
mode properties, when there is a small frequency distance 
between modes with the same azimuthal order m and spher- 
ical harmonic degree 1 differing by 2. In this case each of the 
coupled modes must be represented by a certain superposi- 
tion of spherical harmonics of the two models involved. This 
effect was invoked to explain the nonradial character of pul- 
sation in /3 Cep stars but it could be a possible motivation 
for the worse match of the three last frequencies, either with 
the closest frequencies with 1=2, or with those with 1=0. It 



Table 2. Observed frequencies of Star 278. In the third and fourth 
column we have reported the 1 and n identification of the best 
fit model obtained searching inside the observational box. With 
asterisks are labelled those frequencies which are not matched. 
For this frequencies we report 1 and n of the nearest frequency 
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is important to stress that this is only a suggestion and that 
a definitive statement will be possible only when better data 
will be available. 



5 EXAMPLE-2: STAR 278 OF NGC 6530 

Now we analyze Star 278 in NGC 6530. Among the nine 
frequencies emerged from the Fourier analysis of the data 
performed by ZW07, G07 were able to fit aU with 1=0,1,2 
p modes, except the lowest one. All their best fit models lie 
at a much lower T^b (1800K below the mean temperature 
given by ZW06). The authors support this interpretation 
considering that this star, together with some other stars in 
the cluster, might have been dereddened in excess in ZW06, 
when their color temperatures were derived. To support this 
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hypothesis, they notice that the amount of gas near these 
stars correlates linearly with the difference between model 
and observed surface temperature. 

This Toff is really very low, if comp a red to the red edge of the 
IS calculated by iMarconi fc Pallal (|l998l ) (~ 600K see Fig. 
[6} . As pointed out in Section 3 the cool edge of the radial IS 
also applies to non radial oscillations, and it only depends 
on the efficiency of convective treatment, being cooler for a 
smaller efficien cy of convection. In part icular, the red edge of 
the radial IS of iMarconi fc Pallal l|l998l ) was obtained with a 
value of 1/Hp=1.5 and its effective temperature decreases by 
~ 200 K for 1/Hp = l remaining still hotter than the best fit 
of G07. The efficiency of convection in pre-MS can be tested 
for the convective Hayashi phase of the low mass stars by 
using the location in the HR diagram of pre-MS binaries. 
The binaries however do not provide a calibration for the 
stars located on the 'radiative' phase of the tracks, which 
is practically independent on the convection model. In fact, 
the analysis of dynamical PMS masses shows a good agree- 
ment with all stellar models for M^1.2M0 (Hillenbrand & 
White 2004) for which there are no stars on the convective 
tracks in the examined sample. We know that at the solar 
Tefi (5780K) the l/Hp ratio must be quite large for models 
to fit the solar radius. We also have further hints for mod- 
els at larger Toff, where the convective layer becomes very 
thin. Non local models of convection (Kupka & Montgomery 
2002) for A stars between 7200 and 8500 K can not be eas- 
ily compared to MLT results, but comparable values of the 
maximum convective flux in the H-convection zone are ob- 
tained with a=l at Tofj ~ 7100K (log Tofr=3.85) and with 
a=OA at Tcff ~ 8000K (log Tofr=3.90). Thus non local mod- 
els shows that the convection efficiency is much smaller than 
in standard MLT models at the blue boundary of the IS, but 
rapidally increase while convection deepens. It is clear that a 
further exam of the IS boundaries is needed. For the present 
test we prefer to search in the observational box allowing 
only temperature s hott er than the IS red edge obtained by 
IMarconi fc Pallal (|l998l ). 

We have studied this star using our grid of non rotating 
models since no evidence of relevant rotation is observed in 
the spectrum. This does not mean that the star does not 
rotate, but it could mean that the star rotates so slow that 
the effect on the frequencies is smaller than the observational 
uncertainties. Using the luminosity and temperature given 
by ZW06 and the same size of the observational box used 
for Star 85 we are able to reproduce only seven of the nine 
frequencies. It is important to say that among the models 
able to reproduce seven frequencies only one model in the 
grid has <1: the position of this best fit model in the 
HR diagram is reported in Fig. [5] and the mode identifi- 
cation in Table 2. The only non matched frequencies are 
those two which lie very close to each other. As we noticed 
discussing the results for Star 85 also in this case the the- 
oretical frequencies which are close to these two observed 
ones have / that differs by 2 (see rows marked with asterisks 
in Table 2) and they may emerge in the observed spectra 
as coupled modes due to the rotation, even if slow enough 
that w e don't see it as split of frequencies w ith 1 7^ 0, as sug- 
gested [Daszxns^^Daszkiewicz^^e^^ ||2002| V But a definitive 
conclusion can be done only when new data will be avail- 
able. In this direction also goes the observation that the two 
frequencies that we are not able to match with the theo- 



retical ones are so close, suggesting the presence of a third 
one, that could be a part of a triplet due to the rotation. 
Unfortunately, is impossible to find this frequency from the 
available data. 



6 CONCLUSIONS 

In this paper we have discussed the potentialities of a 
method built up for the interpretation of the observed 
multi-frequency spectra of pulsating stars, giving us the 
possibility of constraining its mass and age. 
The method relies on the direct comparison of the observed 
oscillation spectrum of a given star with a database of 
theoretical ones. How well a given theoretical spectrum 
matches with the observed one is measured by ■ Each 
frequency set in the database was computed using the 
Lose adiabatic oscillation code to derive the pulsation 
frequencies of the evolutionary models computed with the 
ATON code. The characteristics of the grid were reported in 
Section 3 and it is fine enough to determine a best fit model 
to the oscillation data as constrained solely by observed 
frequencies. Till now, the method has been applied only 
to the interpretation of the PMS pulsating stars, but the 
intention is to extend it also to later evolutionary phases. 
We have included the rotation of the star not only con- 
sidering the removal of the m-degeneracy, but also in 
the calculation of the evolutionary models from which we 
compute the pulsational frequencies. During the evolution, 
the angular momentum distribution assumes rigid rotation 
in the convective layers, and conservation of the momen- 
tum shell b y shell , in the radiative core, as suggested by 
IWolff et all (|2004 ). We have shown that, for the some 
physical inputs, the variations of frequencies with rotation 
are small but significative expecially when better data 
will come, from current space missions for example. In 
particular, our calculations have demonstrated that the 
differences increase when the rotation velocity increases. 
While the frequencies of the rotating stars are smaller 
than those computed for non-rotating models for n>5-6, 
inversely for small n the oscillation frequencies of the 
rotating models are larger then the non rotating ones, up 
to 2-3% for stars that rotate with Voq 80 km/s. 
We have tested the mode-matching method on data of 
two stars of NGC 6530. For the Star 85 we confirm with 
self-consistent rotating models the interpretation of the 
data performed by G07, attributing three close frequencies 
to the mode n=4, 1=1 and m=0, -1-1, -1. We also found that 
the remaining three frequencies are better matched when 
rotating models are used. For the Star 278 we find a different 
fit for the frequencies, corresponding to a model within the 
original error box of the star, and dispute the possibility that 
this star has a TeS much cooler that the red boundary of the 
radial instability strip published by Marconi & Palla (1998). 



ACKNOWLEDGMENTS 

It is a pleasure to thank R. Scufiaire and J. Montalban 
for having provided us the Liege OScillation Code (LOSC) 
and the referee W. W. Weiss for his useful suggestions that 



Matching the frequency spectrum of PMS stars. 



9 



improved the final version of the paper. 
Financial support for this study was provided by MIUR 
under the PRIN project "Asteroseismology: a necessary 
tool for the advancement in the study of stellar structure, 
dynamics and evolution", P.I. L. Paterno. 



REFERENCES 

Bohm-Vitense, E. 1958, Z. Astroph., 46, 108 
Breger, M. 1972, ApJ, 171, 539 

Canuto, V. M., Goldman, I., & Mazzitelh, I. 1996, ApJ, 
473, 570 

Cunha et al. 2007, &ARv, 14, 217 

Daszynska-Daszkiewicz, J., Dziembowski, W. A., Pamyat- 

nykh, A. A., & Goupil, M.-J. 2002, A&A, 392, 151 
Kurtz, D. W., & Marang, F. 1995, MNRAS, 276, 191 
Donati, J.-F., Semel, M., Carter, B. D., Rees, D. E., & 

Cameron, A. C. 1997, MNRAS, 291, 658 
D'Antona, F., Cardini, D., Di Mauro, M. P., Maceroni, C, 

Mazzitelli, I., & Montalban, J. 2005, MNRAS, 363, 847 
Endal, A. S., & Sofia, S. 1976, ApJ, 210, 184 
Ferguson, J. W., Alexander, D. R., AUard, F., et al.2005, 

ApJ, 623, 585 

Grevesse, N., & Noels, A. 1993, in Origin and Evolution 

of the Elements, ed. N. Pratzo, E. Vangioni-Flam, & M. 

Casse (Cambridge Univ. Press), 15 
Guenther, D. B. & Brown, K. I. T. 2004, ApJ, 600, 419 
Guenther, D. B.,Kallinger, T., Zwintz, K., Weiss, W. W., 

& Tanner, J. 2007, ApJ, 671,581(G07) 
Herbst, W., Mundt, R. 2005, ApJ, 633, 967H 
Hillenbrand, L. A. & White, R. J., 2004, ApJ, 604, 741 
Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 
N. R. Landin, P. Ventura, F. D' Antona, L. T. S. Mendes, 

L. P. R. Vaz, 2006, 456, 269 
Lovekin, C. & Deupree, R., 2008, ApJ, in press 

(arXiv0803.2053) 
Marconi, M., & Palla, F. 1998, ApJl, 507, L141 
Marconi, M., Ripepi, V., Palla, F., & Ruoppo, A. 2004, 

CoAst, 145, 61 
Mendes, L. T. S., D'Antona, F., & Mazzitelh, I. 1999, A&A, 

341, 174 

Petit P., Donati J.-F., Cameron A., 2002, MNRAS, 334, 
374 

Ripepi, v.. Marconi, M.. Palla, F., Bernabei, S.. Ruoppo, 
A.. Cusano, F.. Alcal, J. M., 2006a, MmSAI, 77, 317 

Ripepi, v., Bernabei, S., Marconi, M., Palla, F., Arellano 
Ferro, A., Bonanno, A., Ferrara, P., Frasca, A., Jiang, X. 
J., Kim, S.-L. et al. 2006b, A&A, 449, 335 

Ripepi, v., Bernabei, S., Marconi, M., Ruoppo, A., Palla, 
F., Monteiro, M. J. P. F. G., Marques, J. P., Ferrara, P., 
Marinoni, S., Terranegra,L. 2007, A&A, 462, 1023 

Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, ApJ, 
456, 902 

Royer, F., Gerbaldi, M., Faraggiana, R., Gmez, A. E. , 

2002, A&A, 381, 105(a) 
Royer F., Grenier S., Baylac M.-O., Gmez A. E., Zorec J., 

2002, A&A, 393, 897(b) 
Ruoppo, A., Marconi, M., Marques, J. P., Monteiro, M. J. 

P. F. G., Christensen-Dalsgaard, J., Palla, F., Ripepi, V. 

2007 A&A, 466, 261R 



Saumon, D., Chabrier, G., Van Horn, H. M., 1995, ApJs, 
99, 713 

Scufiaire, R., Montalban, J.Theado, S., Bourge, P.-O., 
Miglio, A., Godart, M., Thoul, A., Noels, A. ,2007, 
(farXiv:0712.3474 ) 

Stolzmann, W. & Blocker ,T., 2000, a, 361, 1152 

Treu, T. et al., 2003, ApJ, 591, 53 

Unno, W., Osaki, Y., Ando, H., & Shibahashi, H. 1979, 
Nonradial Oscillations of Stars, (University of Tokyo 
Press) 

Wolff, S. C, Strom, S. E., & Hillenbrand, L. A. 2004, ApJ, 
601, 979 

Ventura, P., Zeppieri, A. D'Antona, A., Mazzitelli, I., 1998, 
a, 334, 953 

Ventura, P., D'Antona,F., & Mazzitelh, I. 2007, Ap&SS, 
420 

Zwintz, K., & Weiss,W. W. 2006, A&A, 457, 237(ZW06) 
Zwintz, K. 2008, ApJ, 673, 1088 

Zwintz, K., Guenther, D. B., & Weiss, W. W. 2007, ApJ, 
655, 342 



